Interference of two photons of different color 
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We consider the interference of two photons with different colors in the context of a Hong-Ou- 
Mandel experiment, in which single photons enter each of the input ports of a beam splitter, and 
exit in the same, albeit undetermined, output port. Such interference is possible if one uses an 
active (energy-non-conserving) beam splitter. We find scenarios in which one "red" and one "blue" 
photon enter the beam splitter, and either two red or two blue photons exit, but never one of each 
color. We show how the precise form of the active beam-splitter transformation determines in what 
way the spectral degrees of freedom of the input photons should be related to each other for perfect 
destructive interference of the different-color components in the output. We discuss two examples 
of active beam splitters: one is a gedanken experiment involving a moving mirror and the other is 
a more realistic example involving four-wave mixing in an optical fiber. 



I. INTRODUCTION 

Bosons, in contrast to fermions, tend to occupy the same state when given the opportunity. This tendency arises 
from the fundamental quantum commutation relations of bosonic operators, and is responsible, for example, for the 
creation of laser light and the Bose-Einstein condensate (BEC). In the case of coherent matter waves, such as those 
in a BEC, stimulated scattering of bosonic atoms enhances the probability for the number of atoms in an already 
occupied atomic mode (in particular, the ground state of their external motion) to increase with time, and to grow 
to macroscopic proportions. 

At the level of a few individual photons (as opposed to stimulated emission in lasers, in which the presence of many 
photons is crucial), the propensity of two photons to bunch together was first observed experimentally by Hong, Ou 
and Mandel (HOM) [1]. Specifically, when two separate photons are incident simultaneously upon a 50/50 beam 
splitter from opposite sides, the two output modes that emerge from the beam splitter each contain either two or zero 
photons, but never one photon each. 

▲ 
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FIG. 1: Hong-Ou-Mandel interference at a beam splitter: Two single-photon wave packets (Pi and P2) impinge on a semi- 
transparent mirror (M) and the output wave packets are detected (D\ and D2). 

Let us illustrate this effect in a simple case, which is illustrated in Fig. 1. (The complete theory is given in 
subsequent sections.) Two input modes (1 and 2) contain one photon each with the same center frequency, temporal 
width and polarization, which arrive at the beam splitter at the same time. The initial two-photon state is denoted 
by 

l*>in = aL4in|vac), (1) 

where ajj^^ and a2j„ are creation operators for input modes 1 and 2, respectively. (These modes need not be monochro- 
matic, but can be wave packets [2, 3].) The relation between the input and output mode operators can be written 
as 

•^lin ~ ''"'^lout ~ P'^LutJ 
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4in = P«lout+™2out> (2) 

where the coefficients r (transmissivity) and p (reflectivity) were assumed to be non- negative, without loss of generality, 
and satisfy the auxiliary equation t"^ + p^ = 1. (Usually, one would give the output operators in terms of the input 
operators, as in the Heisenberg picture, but for our purposes the inverse transformations are more useful.) The input 
and output operators satisfy the commutation relations [ifejOiJ.] = 1 for k = 1 or 2. All other commutators are zero. 
By combining Eqs. (1) and (2), and using the fact that the vacuum is invariant under the action of the beam splitter, 
one finds that the output state is 

|*)out = (r^^ut -p4out)(paiout +™2out)|vac) 

= {rp[{al^tf - (4out)'] + ir' - p')4out4out}|vac). (3) 

In the case of exactly 50% beam splitting, for which t = p = l/V^, Eq. (3) reduces to 

|*)out = (|2)|0)-|0)|2))/x/2, (4) 

where |n) denotes an n-photon Fock state. (All the other modes are still vacuum modes.) In this case the part of the 
state with one photon in each output port is zero. The physical interpretation of this cancelation is that there are 
two paths leading to the same final state component |1)|1): Either both photons transmit through the beam splitter 
or both photons reflect. The probability amplitudes for these two events have opposite signs [r^ and — in Eq. 
(3)], so they cancel. The resulting two-photon bunching is a basic property of nature and is at the heart of quantum 
information processing schemes using linear optics (beam splitters, polarizers, etc.) [4]. 

Notice that this type of HOM interference does not depend on any relative phase between the input states: Unlike 
coherent states or classical light waves (whose interference properties at a beam splitter certainly depend on the relative 
phase of the complex amplitudes), Fock states do not possess any physically relevant phase on which interference could 
depend. 

In the original setting [1], for the HOM effect to work, the photons must be in states that are identical in their 
polarization and spectral degrees of freedom (but their input states differ in their propagation directions). The 
question arises: Must two photons (or other bosons) have identical properties before they come together for this type 
of interference to occur? For example, assuming they have identical polarization states, must two photons initially be 
in identical energy (spectral) states in order to interfere? 

In this paper, we show that the answer is negative: By using an active, frequency-shifting beam splitter, one can 
observe two-photon interference between photons of different color. The active process can be either reflection by a 
moving mirror (for small frequency shifts), or three- wave mixing in an optical crystal [5-8] or four- wave mixing in 
an optical flber (for large frequency shifts) [9-12]. The HOM effect is due to destructive interference in the output 
state, not the input state. We derive conditions on the forms of the input states that are required for high-visibility 
two-photon interference. These phenomena, in addition to demonstrating a fundamental property of light, will have 
applications in quantum information experiments involving photons of different color [13-16]. 



II. PASSIVE BEAM SPLITTERS 



Henceforth, we will use a simplified notation: Instead of writing relations that define the input operators in terms 
of output operators, as in Eq. (2), we will simply write a\ i— >■ Ta\ — pa\ etc., to reflect how the input state changes 
to the output state in the most straightforward way. The operators a\ and appearing on the right-hand side are 
thus meant to be output operators. 

If the beamsplitter is a passive device that preserves the energy and spectrum of each photon separately, then the 
unitary transformation between the input and output channels can be written as 

a\{Lo) I— >■ Ta|(w) — /9a2(w), 

4('^) pa[{u)) + Ta2{oj)- (5) 

Here we use creation operators a\{uj) for monochromatic modes with frequency uj impinging on one input port of 
the beam splitter and a\{u)) for modes impinging on the other input port. (We exploit a quasi one-dimensional 
picture of propagating single-photon wave packets. This picture is valid in the paraxial limit and when the radiation 
is sufficiently narrowband [2, 3].) These operators satisfy the commutation relations [afc(u;), a|.(w')] = 5{ijJ — to') for 
fc = 1 or 2. All other commutators are zero. Suppose now we have an input state containing two photons, one 
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impinging on each input port, described by spectral amplitudes (pi{u)i) and ^2(^2), respectively, which satisfy the 
normaUzation conditions 



^ doJk\(i>k{(^k)f = 1- (6) 
Then the input state is 

|'S')in = j cLi;i(/)i(wi)a|(wi) j duj2<l>2{i^2)al{L02)\vac}. (7) 

This input state is mapped by the beam-splitter transformation onto the output state 

|*)out = j J dujiduj2(l>i{^i)4'2{^2)[rpa\{uji)a\{uj2) - Tpal{ui)al{L02) 

+ r^aj (wi)a2(w2) — p'^a\{uj2)al{uJi)]\va.c). (8) 

The last two terms can be made to cancel each other under the right conditions, and that complete destructive 
interference would correspond to HOM interference. The interference is perfect whenever = and 

(/>l(wi)02(w2) = </'l(tJ2)'/>2(wi) (9) 

for all frequencies wi and UJ2- This condition can be rewritten as 

(/>l(wi)/(/)2(u;i) = ^il(w2)/02('^2)- (10) 

The left-hand and right-hand sides of Eq. (10) are functions of different variables, implying that they cannot actually 
depend on those variables. Hence, 

M^i) = cM^2), (11) 

where C is a constant. Normalization of the quantum states gives |C| = 1, so C is just the phase factor exp{i9) for 
some phase 6. Equivalently, in the time domain, the Fourier transforms of the spectral density functions must satisfy 

Mt) = cMt) (12) 

for destructive interference. 

From these conditions one may get the impression that it is important to have indistinguishable photons at the 

input. But that is a little misleading; what counts is indistinguishability at the output, as we will show. The same 
output state can be reached by two different paths, and it is the interference between the two paths that matters (as 
always in quantum mechanics). 



III. ACTIVE BEAM SPLITTERS 
A. Moving semi-transparent mirror 

Consider a semi-transparent mirror moving to the left in one dimension. In the mirror frame, two channels with 
the same carrier frequency (color), and opposite propagation directions (right and left) impinge on the mirror [Fig. 
2(a)]. The input and output channels are related by Eqs. (5). However, in the laboratory frame [Fig. 2(6)], light 
initally propagating to the right will be blue-shifted upon reflection, whereas light initially propagating to the left will 
be red-shifted upon reflection (and transmitted light does not change frequency). 

It is convenient to label the right-propagating channel "red" and the left-propagating channel "blue." (Of course, 
in order to shift actual red light to actual blue light and vice versa, the speed v would have to be of order c.) The 
moving mirror transforms mode operators in the following way: 

ajj(w) ra|j(w) - (p/Va)aya;/a), 

a^(w) I-)- {p^/a)a^J^{au}) + Ta^g{u)). (13) 
The constant a is determined by the Doppler effect, which acts twice: 



a=(l-/3)/(l + /3)<l. 



(14) 
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FIG. 2: Interaction of two waves (r and /) at a semi-transparent mirror (M). In the mirror frame, the waves have the same 
frequency and opposite directions of propagation (right and left). In the laboratory frame, the waves have different frequencies 
and directions. 



where the normalized mirror speed is /? = v/c. The factors of a in Eqs. (13), which are absent from Eqs. (5), are 
required to conserve the number of photons, because red modes in the frequency interval dojR are coupled to blue 
modes in the interval dujB = dun/a. Transformation (13) is unitary, preserving the commutation relations for the 
transformed operators. For example, 

[a_B(aj), a}g{Lo')] H- [py/aaji{auj) + rasii^), p^/aa]^{auj') + ra^(aj')] = S{oj — uj'), 

[aB{uj), a^(Q!w')] [py/aaRiauj) + rasiuj), Ta^j^{auj') - {p/ ^/a)a\{uj')] = 0, (15) 

where we used i5[Q;(a; — w')] = (5(w — uj')/a. Although this transformation preserves the number of photons, it does 
not preserve their energy, thus earning the right to be called active. 

Now suppose we have one photon in each input mode, described by wave packets centered around red and blue 
frequencies, as the input state: 

\'i)-in = j diuii(l)B.(u)R)a\{LUii) j diUB4>B{'^B)a^B^WB)\vac), (16) 

where the spectral amplitudes satisfy the normalization conditions (6). This state is transformed by the moving 
mirror into 

|*)out = j J dujRdujB<j)R{ujB.)4'Bii^B)[Ta^B.iujR) - ip/Va)a^Bi^ii/'^)] 

x[/3Vaaij(awB) + Ta^(a;s)]|vac). (17) 

The interesting effect is that there can be destructive interference in the output terms where one red and one blue 
photon emerge. There are two amplitudes for reaching that final state with opposite signs; the (necessary and 
sufhcient) condition for exact destructive interference for all frequencies ujr and lob is 

t^(I)r{ll>r)(I)b{(^b) ^ p^(t)R(auB)4'B{'^R/a)- (18) 
This condition implies that = p^ and 

4'R{<^R)/4'B{^^R/ct) = (f>R{aUJB)/<pB{i^B) (19) 

for all frequencies uir and ujb, which in turn implies that 

(t^Ri^R) = C(f>B{ujR/a), 

M^b) = C-^Ma^B), (20) 

with C a constant independent of ujr and lur- These two conditions are really the same. Normalization gives 
|Cp = 1/a, so we can write C — e'Kp{i0) / ^/a, with 6 some constant phase. In the time domain this condition gives 



4>Bit) -.^ J dujB4>Bii^B)exp{iu}Bt) ^ C J d{ujR/a)4>R{u}R)exp{iu}Rt/a), (21) 

which means that 

V^exp(i6')0B(t) = ^Rit/a). (22) 

If the input wavepackets satisfy these relations, then the output state of the moving mirror will contain either two 
red photons moving to the right, or two blue photons moving to the left, but never one of each color. 
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B. Frequency translation in fiber, simple version 



Now consider four-wave mixing in a fiber. Among the many processes that may occur, we are interested in frequency 
translation (Bragg scattering), in which two strong, classical pump waves (p and q) convert "signal" (s) photons into 
"idler" (r) photons, or vice versa [9-12]. The frequency-matching condition for this process, which is illustrated in 
Fig. 3, is 



(23) 



It follows from Ec^. (23) that the difference between the signal and idler frequencies, — Wr — equals the difference 
between the pump frequencies, ujg — ujp. It is convenient to refer to the idler as the "red" sideband and the signal 
as the "blue" sideband. (Once again, we use these terms just to indicate two different colors.) If the pumps are 
monochromatic (CW), then each red frequency is coupled to one blue frequency. If the red frequencies all lie within 
the phase-matched bandwidth of the convertor, the red-to-blue conversion efficiency is frequency independent. If the 
pump powers and fiber length are chosen judiciously, the probability for frequency translation to occur for a given 
photon is 50%, with the remaining 50% probability assigned to the photon's frequency staying the same. 
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FIG. 3: Two classical pumps (p and q) can drive a four-wave mixing process in which a single signal (s) photon is frequency- 
converted to a single idler (r) photon, or vice versa. 



The transformation between the red and blue modes is [9, 10] 

a^(wfl) i-^ Ta\^{uR) - pa^g{uR + rt), 

agiujB) pa^ws - il) + ra^(wB). (24) 

(An equivalent transformation governs frequency up-conversion in a crystal [5].) This transformation is unitary. By 
following the steps between Eqs. (16) and (18), one finds that the HOM interference is perfect when 

T^(l)ii{u>ji)(f>B{LL>B) = p^(I)r{ujb - fl)4)Bi(^R + (25) 
Condition (25) implies that = and 

4>r{^r) = C4)b{c^r + ^), 

M^b) - C^^M^B-n). (26) 

with C a constant, independent of ujr and lob- Once again, these two conditions are the same. Normalization gives 
\C\ = 1, and hence one can write C = exp(i0). In the time domain, 

eM^0-i^t)4>B{t) = 4>R{t), (27) 

which means that the temporal shapes should be identical, with an overall frequency difference equal to f2. Equation 
(27) provides necessary and sufficient conditions for frequency translation to display HOM interference in this simple 
case: Two input photons, one red and one blue, traversing a nonlinear fiber with the appropriate amount of frequency 
translation, lead to an output state which contains either two red photons or two blue photons, but never one of each 
color. 



C. Frequency translation in fiber, general version 



If one of the pumps is pulsed, or both pumps are pulsed, each frequency component of the signal is coupled to many 
frequency components of the idler and vice versa [24]. In either case, the transformation over all frequencies is of the 
form 

a\uj)^ I dLo'G{uj,uj')a\Lo'), (28) 
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where the kernel G satisfies the unitarity condition 

du'G{io,uj')G*{co",uj') = 5{oo - cj"). (29) 



Let us split the frequency interval into two non-overlapping parts, "red" and "blue," and let us correspondingly write 
a{u)) = aii{u)ii) whenever a; is a "red" frequency, and a{uj) = asioJE) whenever a; is a "blue" frequency. Then we can 
rewrite the kernel in the block-matrix form 



Grr{u:r, uj'j^) Grb {'^r- '^''b) 
GbrX^Bt^'r) Gbb{^b,^''b) 



(30) 



where the constituent (Green) kernels Grr, Grb, Gbr and Grb describe red-red transmission, blue- red conversion, 
red-blue conversion and blue-blue transmission, respectively. The Green kernels satisfy the unitarity conditions 

j duiRGRR{UR,0jR)G*RR{LOR,UJR) + j duiBGRB{ujR,L0'B)G*RB{^R,^'B) = S{lor-ujr), 

J rfa;^GBK(wB,a;^)GJj^(Wfl,a;^) + j rfa;^GBB(a;B, wJj)GJjs(w^,a;^) = 0, (31) 

and similar conditions in which R -(r^ B. In terms of these kernels, the transformations between the red and blue 
modes are 

Oiiii^R) I-*- j dx^RGRR{ujR,u}R)a^ji{uj'R) + J dLj'BGRB{i^R,u;'B)aBi^B), 

05(^5) I-*- j dwijGBii(wB,Wfl)a^(Wfl) + j dio'BGBB{oJB,i^B)aB{^B)- (32) 
The interference of one red and one blue photon in the final state is perfectly destructive when 



// 



duiRduJB 4'R (Wfl ) (f>B (WB ) [Grr {ur ,uj'i{)Gbb (wb , Wb ) 

-I- GRBii^R;L-j'B)GBR{(^B,(^'R)] = (33) 



for all frequencies and ui'q. By rearranging the terms in Eq. (33), so that functions of and ui'g appear on 
difl[erent sides, one obtains the interference conditions 

j du)R(l)R{uiR)GRR{uiR,Ul'ji) = ^ j duiB4'B{0JB)GBR{0JB,0jR), 

j diOB(j>B{<^B)GBB{i^B,i^'B) = j '^R^r{'^r)^Rb{^R,^'b)^ (34) 

where C is a constant, independent of frequency. The first condition states that the spectral amplitude of red output 
photons arising from red input photons be equal (up to a constant) to the spectral amplitude of red output photons 
arising from blue input photons. The second condition is similar, but for blue output photons. 

In general, though, these conditions will not be fulfilled with given pumps in a given fiber. One needs to design the 
four-wave mixing process in such a way that with 50% probability a red photon is cither converted to a blue photon 
or remains a red photon with a possibly altered spectral wave function. The opposite transformations must hold for 
the incoming blue photon. The fiber dispersion, pump pulses, and the input red and blue wave packets must all be 
designed properly to ensure that the output red and blue wave packets are identical for either input channel. (This 
is work in progress.) 

More precisely, if one uses the unitarity conditions (31) to rewrite and combine Eqs. (34), one obtains the interference 
conditions 

dijJBK{ljOR,bJB)<f)B{l^B), 
(t>B{u)B) = j dWRK*{ujR,UJB)(t)R{i^R), (35) 

where the HOM kernel K{ur,cob) is 

K{ujr,ujb)= j du'[G*RR{ujR,uj')GBR{uJB,i^') - G*rb{^R^^')Gbb{u)b,co')]. (36) 
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Recall that every hermitian kernel can be decomposed in terms of its eigenvalues, which are real, and a single set 
of eigenfunctions, which are orthonormal. In a similar way, every complex kernel has the Schmidt (singular-value) 
decomposition [17, 18] 

K{lOr, UJb) = J2 (^nRnMBliuJB), (37) 
n 

where the coefficients (singular values) (t„ arc real and non-negative, and the mode functions {_R„(u;/f)} and {_B„(wb)} 
are self-orthonormal. Subroutines are available, which determine the mode functions numerically. By substituting 
decomposition (37) into Eqs. (35), one finds that the interference conditions can be fulfilled only if there is at least 
one unit coefficient among the Schmidt coefficients. If there is exactly one unit coefficient, say <Ji = 1, then one 
must have cf)R{ujji) = exp{i9)Ri{ujii) and (j)B{i^B) = Bi{ujb) in order for complete interference to occur. If there are 
multiple imit coefficients, one can take superpositions of the corresponding modes. If the kernel possesses only unit 
coefficients, then one can choose any red spectral amplitude 4'r{(jOr) and find a corresponding blue spectral amplitude. 
The simple cases of the preceding subsections are, in fact, such degenerate cases, whenever = p^- In the general 
case, which corresponds to pulsed pumps, one expects that only certain red and blue wavepackets, of finite bandwidth 
and duration, lead to perfect HOM interference. In all of these cases, C must again be the phase factor exp(z0). Thus, 
the fiber and pump requirements for perfect HOM interference may be stated in a different and more succinct way: 
The HOM kernel must have at least one unit singular value. 

One can analyze the general case in an alternative way, by giving the Schmidt decompositions of the Green kernels 
directly, rather than the HOM kernel K. By using the results of [19] , which are reviewed in the Appendix of this 
paper, one finds that 

Grr{uJr,Jr) = ^TnVn{y>R)vl{u'R), 

n 

Grb{^^R,^^'b) = -^PnVn{y)R)wl{uj'B), 
n 

Gbr{ojb,oj'r) = ^p„W„(wbX(w^), 

n 

Gbb{u1B,U]'b) = ^TnWn{uiB)w^{uj'B), (38) 
n 

where the non- negative coefficients t„ and /)„ satisfy the auxiliary equations + Pn = li a-nd the sets of mode 
functions {Vn}, {vn}, {Wn}, and {«;„} are self-orthonormal. By combining Eqs. (36) and (38), one obtains the 
alternative decomposition 

K{ur,UJb) = 2Y,PnrnV:{iJR)Wn{uJB). (39) 

n 

It follows from Eqs. (37) and (39) that cr„ = 2p„T„, i?„ = V* and B„ = W*. Notice that a„ < 1, so the HOM 
condition (7„ = 1 represents a special (ideal) case. 

The coefficients r„ and pn can be interpreted as generalized beam-splitter coefficients, with r„ playing the role of a 
transmission coefficient (not changing the color) and pn the corresponding reflection coefficient (changing the color), 
for each mode n. Specifically, the Schmidt-mode expansions 

al{(w)|in = ^atl4(w), 

n 

n 

ai{(w)|out = ^civn{u:), 

n 

4 ('^) lout = ^diwn{uj), (40) 

n 

allow one to rewrite the continuous input and output operators for red and blue light in terms of the discrete operators 
{an,bn,Cn,dn} [2]. In terms of these specially constructed discrete modes, the output operators are related to the 
input operators in the simple, pair-wise manner 

bi = Pnci + Tndl (41) 
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Equations (41) are equivalent to the two-mode equations (2). As expected, a unit ij„ corresponds to a 50/50 beam- 
splitter transformation with Tn = Pn = 

Moreover, description (38) allows us to quantify the effect of imperfections in frequency translation. That is, suppose 
that the fiber process we actually designed does not have any coefficient t„ equal to 1/a/2- Then the best we can do 
is to use the mode that corresponds to the transmission coefficient that is closest to l/-\/2, say Ti « l/-\/2. In fact, 
if we choose in that case ^r{u)r) = Vi{ujr) and ^s(wb) = Wi{ujb), then the amplitude of the component of the 
output state containing one red and one blue photon is 

(vac|afl(wfl)as(a;B)|*)out = (rf - pl)V{{uR)W*{iJB). (42) 

This equation generalizes the corresponding term in Eq. (3). The deviation from perfect two-photon interference can 
be quantified as the probability to find one red and one blue photon (at any red-blue frequency pair) in the output, 
and (42) shows that this probability is Prb = {tI — pf)^- 



IV. CONCLUSIONS 



In this paper, we showed that one can observe^ interference between two photons of different color. These photons 
arise independently from separate sources before they come together and interfere. The important point is that both 
photons end up in the same final spectral state (or mode), even though they began in two spectrally distinct states. 
It is not necessary for them to be identical initially, except in the obvious sense that they are both photons, which are 
fundamentally identical particles (if one takes a particle view of photons). For example, if one "red" and one "blue" 
photon enter the active interferometer, then either two red or two blue photons exit, but never one of each color. 

Observing such bosonic bunching of photons, whose colors are different initially, requires reversible and symmetric 
inelastic scattering processes. Examples of such processes include three-wave mixing (wavelength down-conversion) 
in a crystal [6-8] and four- wave mixing (Bragg scattering) in a fiber [11, 12]. Both types of process can shift photon 
wavelengths by several hundred nanometers, with high photon-conversion efficiencies. Such two-photon interference 
should find applications in quantum information science, in schemes using entanglement between optical modes of 
different color. 
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Appendix A: Decomposition of a Unitary Transformation 



In this appendix, it is shown that a 2n-mode unitary transformation is equivalent to n 2-mode beam-splitter-like 
transformations. As an illustrative example, consider the effects of a rotator (circular phase-shifter) on a monochro- 
matic wave [20, 21]. Let and A- be the components of the 2x1 amplitude vector A relative to the circularly- 
polarized basis vectors = [l,i]*/2^/^ and E- = [1, — i]*/2^/^, respectively, and let 6 be one-half of the differential 
phase shift. Then the input and output vectors are related by matrix equation 



A{e) = M{e)A{o), 



(Al) 



where the 2x2 transfer matrix 



M((9) 



e'^ 
e-'' 



(A2) 



is diagonal, with complex-conjugate eigenvalues. Now let A^ and Ay be the amplitude components relative to the 
linearly-polarized basis vectors Ex = [1,0]* and Ey = [0, 1]*, respectively. Then the associated transfer matrix 



M{0) 



cos 6 sin 
— sin 6 cos 6 



(A3) 
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What appears as phase shifts relative to one basis appears as rotation relative to another. The second transfer matrix 
also describes the effects of a beam splitter on two input waves (t = cos 9 and p = sin 9) . 
Now consider the 2n x 2n unitary matrix 



M 



A B 
C D 



(A4) 



where the n x n matrices A, B, C and D describe signal-signal transmission, idler-signal conversion, signal-idler 
conversion and idler-idler transmission, respectively. Our derivation of the beam-splitter decomposition of a unitary 
transformation is similar to the derivation of the Bloch-Messiah decomposition of a squeezing transformation [22-24]. 
The unitarity conditions are 



/ 
/ 



(AAt + BSt) {AC'^+BD^) 
{A^A + C^C) {A^B + C'^D) 



The complex matrix A has the singular-value decomposition (SVD) 

A = U^D^Vl 



(A5) 



(A6) 



where Ua and Va are unitary, and Da = diag(aj). In the standard SVD, the aj are non-negative. However, one 
can also allow them to be negative. Similar SVDs exist for the other transmission matrices, in which Df, = diag(/3j), 
Dc = diag(7j) and = diag((5j). The first diagonal condition in Eq. (A5) implies that 



UaDlul + U.Dlul = J, 



(A7) 



from which it follows that Ua = Ub and + ~ 1. The other diagonal conditions imply that Uc = Ud, Va — Vc, 



(5j and /3| 



7? The off-diagonal conditions all imply that 5j = aj and = —fij. Let Eij be the jth 



column of Ua, E2j be the jth column of Ud, Fij be the jth column of Vc and F2j be the jth column of Vb (so Eij is 
an eigenvector of AA'^ , E2j is an eigenvector of DD^, Fij is an eigenvector of C^C and F2j is an eigenvector of B^B). 
Then these vectors form self-orthonormal sets and 



-E^,I3,fI Eo.a.Fl 



(A8) 



Equation (A8) is the beam-splitter decomposition [19]. Fij and F2j are the input eigenvectors, whereas Eij and E2j 
are the output eigenvectors. It follows from Eq. (A8) that 



Fi,ajE\. -Fij^jEl- 
F2jpjEl 



F2jajElj 



(A9) 



Thus, if the forward transformation is determined by the parameters (aj, (3j, — /3j, aj), the backward transformation 
is determined by the parameters {aj,^j3j, I3j,aj), just like the two-mode transformation (A3). Furthermore, the 
roles of the input and output eigenvectors are reversed. Equations (AS) and (A9) remain valid if one replaces the 
aforementioned real parameters by the complex parameters (aj, — /5*, a*) and [a*. —0j, 0* ,aj) in the forward and 
backward transformations, respectively, where jajj^ + = 1- 

The Green kernels of the text are related to the constituent matrices of this appendix by the relations Grr = A*, 
Grb = C*, Gbr = -B* and Grr = -D*, from which it follows that the HOM kernel K = A^B - C^D. By combining this 
result with decomposition (A8), one finds that 



K = 2j2a*PjFijF^.. 



(AlO) 



Equation (AlO) is consistent with Eq. (39), which defines the mode functions of the text {(l)r and (pb) indirectly. These 
mode functions are defined directly by the interference conditions (35). In the matrix notation of this appendix, they 
are the eigenvectors of KK^ and K^K, respectively. By combining the formula for K with conditions (A5), one finds 
that 



AA'^AC^C, 
AB^fBD^D. 



(All) 
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a'' A and C^^C have the eigenvalues and |/3jP, respectively, and the (common) eigenvectors Fij, whereas B^B 

and D^D have the eigenvalues |/3jp and la^p, respectively, and the (common) eigenvectors -Fjj- Hence, the mode 
functions (f>rj — Fij and (pbj = F2j. The associated singular values = 2\a jl3j\ are the square roots of the (common) 
eigenvalues of KK^ and K^K. These results are consistent with Eq. (37). 
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